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ABSTRACT 



O 

o 

' Context. Some intra-day variable, compact extra-galactic radio sources show brightness temperatures severely exceeding 10'- K, the limit set 
, by catastrophic inverse-Compton (IC) cooling in sources of incoherent synchrotron radiation. The violation of the IC limit, actually possible 
IL* under non- stationary conditions, would lead to IC avalanches in the soft-y-ray energy band during transient periods. 
. ^ Aims. For the first time, broadband signatures of possible IC catastrophes were searched for in a prototypical source, S5 0716+71. 

Methods. A multifrequency observing campaign targetting S5 0716+71 was carried out during November 06-20, 2003. The observations, 
5_j organized under the framework of the European Network for the Investigation of Galactic nuclei through Multifrequency Analysis (ENIGMA) 
together with a campaign by the Whole Earth Blazar Telescope (WEBT), involved a pointing by the soft-y-ray satellite INTEGRAL, optical, 
near-infrared, sub-millimeter, millimeter, radio, as well as Very Long Baseline Array (VLB A) monitoring. 

Results. S5 0716+71 was very bright at radio frequencies and in a rather faint optical state (i? = 14.17- 13.64) during the INTEGRAL pointing; 
significant inter-day and low intra-day variability was recorded in the radio regime, while typical fast variability features were observed in the 
optical band. No obvious correlation was found between the radio and optical emission. The source was not detected by INTEGRAL, neither 
by the X-ray monitor JEM-X nor by the y-ray imager ISGRI, but upper limits to the source emission in the 3-200 keV energy band were 
estimated. A brightness temperature > 2.1 x lO''* K (violating the IC limit) was inferred from the variability observed in the radio regime, 
but no corresponding signatures of IC avalanches were recorded at higher energies. 

Conclusions. In the most plausible scenario of negligible contribution of the interstellar scintillation to the observed radio variability, the 
absence of the signatures of IC catastrophes provides either a lower limit (5 S 8 to the Doppler factor affecting the radio emission or strong 
constraints for modelling of the Compton-catastrophe scenario in S5 0716+71. 

Key words, galaxies: active - galaxies: BL Lacertae objects: general - galaxies: BL Lacertae objects: individual: S5 0716+71 - galaxies: 
quasars: general - gamma-rays: observations - radiation mechanisms: non-thermal 
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* Partially based on observations obtained with INTEGRAL, an 
ESA project with instruments and science data centre funded by 
ESA member states (especially the PI countries: Denmark, France, 
Germany, Italy, Switzerland, Spain, Czech Republic and Poland), and 
with the participation of Russia and the USA. Partially based on obser- 
vations by the Whole Earth Blazar Telescope (WEBT); for questions 



regarding the availability of the data from the WEBT campaign, please 
contact M. Villata (villataOto . astro . it). Partially based on ob- 
servations collected at the German-Spanish Astronomical Center, 
Calar Alto, operated by the Max-Planck-Institut fUr Astronomie, 
Heidelberg, jointly with the Spanish National Commission for 
Astronomy. Partially based on observations with the 100-m tele- 
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1. Introduction 

The phenomenon of intra-day variability (IDV; Wagner & 
Witzel [T995t has been a long-standing problem since its disco- 
very in the late 1960s (see e.g. Racine et al. 119701 Witzel 
et al. \T^EE[ Heeschen et al. .1987^ . The occurrence of IDV 
appears to be more common in flat-spectrum extragalactic 
sources dominated by a very compact core in Very Long 
Baseline Interferometry (VLBI) maps (Quirrenbach et al. 1992. 
I2000> . With recent Very Long Baseline Array (VLBA) mea- 
surements, Kovalev et al. (2005 1 showed that IDV sources typi- 
cally exhibit a higher compactness and core-dominance on sub- 
milliarcsecond scales than non-IDV ones; they also found that 
a higher amplitude of intraday variations characterizes sources 
with a higher flux density in an unresolved VLBA component, 
and that the most variable sources tend to have the most com- 
pact structure. 

Rapid variations in flux may be caused by mechanisms in- 
trinsic to the source. In this case, causality arguments would 
imply that the variability originates from very compact re- 
gions of the AGN, thus characterized by high photon densi- 
ties and brightness temperatures. In sources of incoherent syn- 
chrotron radiation, very high photon densities would lead to 
catastrophic cooling via inverse-Compton (IC) scattering of 
the synchrotron radiation by the high-energy electrons, with 
a production of high-energy radiation much higher than ob- 
served (Hoyle et al.i966l. Keflermann & Pauliny-Toth (.1969J 
showed that the onset of catastrophic radiation losses, occur- 
ring when the photon energy density in the emission region 
exceeds the energy density of the magnetic field, limits the 
maximum observed brightness temperature to the so-called TC 
Umit' Tb ^ 10'2 K (see also Readhead 1994"). They also no- 
ticed, however, that this limit may be significantly exceeded 
under non-stationary conditions. Evidence of severe violations 
of the IC limit was reported for many compact IDV radio 
sources (see e.g. Quirrenbach et al. .1989. Wagner & Witzel 
119951 Kedziora-Chudczer et al. I1997> . although no informa- 
tion about corresponding IC-scattered emission is available for 
any of them. 

Alternatively, propagation effects, like variations of the ab- 
sorption along the line of sight, deflection of the light in the 
potential well of foreground stars (microlensing; Chang & 
Refsdal ll 9791 1. and, in the radio regime, interstellar scintillation 
(ISS; Rickett 1990 Rickett et al. 1995 1, are possible extrinsic- 
source mechanisms which have been invoked to explain IDV. 
On the one hand, microlensing was shown to be an unlikely 
explanation of the rapid variability in the best-studied IDV 
sources, due to the frequency of the flaring activity, the sta- 
tistical asymmetry of the light curves, the short time scales of 
variability, and the non-zero lag measured between optical and 
radio variations (see Wagner ll 9921 and references therein). On 
the other hand, ISS is not expected to play a major role at mil- 
limeter and sub-millimeter wavelengths, the transition between 
weak and strong scintillation regimes occurring in the centime- 
ter domain (Rickett et al. ll995> . 



scope of the MPIfR (Max-Planck-Institut fiir Radioastronomie) in 
Effelsberg, Germany. 



In order to reconcile intrinsic variations with the theoreti- 
cal limit, different explanations, like beaming of the emission 
due to bulk relativistic motion (Rees U967> . coherent radia- 
tion mechanisms (Baker et al. 1988 Benford|T592 Lesch & 
Pohl 1992 1, propagation of a shock in an underlying, stable 
relativistic jet (Qian et al. I1991> . were proposed. Neither the 
broad-band nature of the spectral energy distribution (SED) 
observed in IDV sources (see e.g. Quirrenbach et al. 119891 
Wagner et al. ll993l|l996> nor the VLBI observations (see e.g. 
Gabuzda l2000> strongly support any of them, leaving the alter- 
native explanation that the IC limit might actually be violated 
and inverse-Compton catastrophes occur during transient pe- 
riods (Kellermann & Pauliny-Toth 1969, Slysh 1992 1. In this 
paper, this hypothesis is explored for a prototypical source, S5 
0716-H7I. 

This source is one of the brightest and best-studied BL 
Lacertae objects in the sky. It was one of the prime targets 
for investigating the mechanism responsible for IDV, and the 
first source in which simultaneous variations in the radio and 
optical bands, indicating a possible intrinsic origin of the ob- 
served IDV, were reported (Wagner et al.'1990' QuiiTenbach et 
al. ll991> . Moreover, it exhibited IDV during all the past optical 
studies and almost all the radio campaigns caiTied out during 
the last two decades (Heeschen et al. 119871 Wagner et al. »l990l 
>1996 Heidt & Wagner 1996 Ghiselhni et al. nWTl Sagar et 
al. 1999 Quirrenbach et al. 2000 Viflata et al.|2000 Nesci et 
al. .2002. Kraus et al. 2003 Raiteri et al. l2003t . The IDV duty 
cycle of the source (the fraction of time in which the object is 
variable) derived from these studies is ~ 90%. 

Deep maps of the source obtained with the Very Large 
AiTay (VLA) show a core-halo structure on the arcsecond 
scale. VLBI observations over more than 20 years at centimeter 
wavelengths reveal a very compact source, with evidence of a 
core-dominated jet structure extending several tens of milliarc- 
seconds to the North (Eckart et al. 119861 |l987 Witzel et al. 
il988; Polatidis et al.fT995l Jorstad et al. lM)li) . The milHarc- 
second jet is misaligned with respect to the VLA jet by ~ 75° 
(Britzen et al. 2005 see also Eckart et al. ll987t . 

Controversial scenarios, involving a wide range of proper 
motions (0.05-1 .2 mas/year), were proposed for the kinematics 
of the S5 0716+71 jet components, the more recent ones poin- 
ting towards apparent velocities which are atypically fast for 
BL Lac objects (see Bach et al. 120051 and references therein). 

The redshift of the source is still unknown, although the 
starlike appearance and the absence of any signature of a 
host galaxy in deep images had set a lower limit of z > 0.3 
(Quirrenbach et al. 1991; Stickel et al. 119931 Wagner et al. 
1996 ), which will be used throughout this paper. More recently, 
Sbarufatti et al. J2005> suggested a higher lower limit of z > 
0.52. The exact brightness temperature of the source can there- 
fore not be determined: lower limits up to T\, ~ lO'^ K were 
inferred from radio IDV at 5 GHz (Quirrenbach et al. I199U 
Wagner et al. 1996), whereas a limit of rb,;=o - 1.85 x 10'^ 
K was estimated from the constraints on the core size derived 
through interferometric measurements performed in August 
2003 with the Very Long Baseline Array (VLBA) at 15 GHz 
(Kovalev et al. 2005 1. This source is hence an ideal target for 
the investigation of radiative signatures of IC catastrophes. 
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S5 0716+71 has been detected at GeV energies with a steep 
y-ray spectrum (Hartman et al. 1999), but the soft-y-ray part 
of its spectral energy distribution is poorly known: upper limits 
to the emission were provided by OSSE (McNaron-Brown et 
al.[l995) and COMPTEL (Schonfelder et al. l2000t : a recent 
reanalysis of the COMPTEL data (Collmar l2006b yielded a 
source detection in the 3-10 MeV energy range. 

Synchrotron-self-Compton (SSC) emission models, which 
can reproduce the GeV y-ray emission of the source (see e.g. 
Ghiselhni et al. I1997> . predict the SED peak to occur in the 
MeV-GeV domain. The MeV detection reported by Collmar 
J200 6^ is indeed consistent with this scenario. The above SSC 
modelling involves a wide range of Lorentz factors of the ra- 
diating particles, up to y ~ 10^, in order to explain the high- 
energy component of the spectrum. 

In a Compton catastrophe, the main role is expected to be 
played by electrons responsible for the bulk of the radio emis- 
sion in the GHz regime, where the violations of the IC limit 
have been observed. Assuming an observed peak frequency 
of the radio spectrum of a few GHz, and a magnetic field 
B S 10^^/(5 fiG (6 being the Doppler factor of the emitting 
region), the above electrons would be characterized by Lorentz 
factors not greater than a few times 10^. First-order IC scat- 
tering by these particles would generate bursts of radiation at 
~ 10'"*- 10'^ Hz, which might blend with the synchrotron radi- 
ation produced by the high-energy tail of the same distribution 
of particles. The dominant loss-term would be second-order IC 
scattering (see e.g. Kellermann & Paulinv-Toth 1 1 9691 Bloom 
& Marscher 1996 1, also responsible for the onset of the ca- 
tastrophe, and would boost photons into the 10'** - 10^" Hz fre- 
quency range. Flares of IC-scattered radiation should therefore 
be observed in this energy range whenever the Compton limit 
is violated. The efficient cooling associated with a Compton ca- 
tastrophe would then rapidly restore the actual brightness tem- 
perature of the source. 

The advent of the soft-y-ray INTEGRAL satellite offers an 
unprecedented chance of investigating this effect. 

A multifrequency campaign involving an INTEGRAL 
pointing of the source and simultaneous radio, millimeter, sub- 
millimeter, near-infrared, and optical photometric monitoring, 
as well as VLBA observations, was organized and carried out 
during November 2003. For the first time, constraints on the 
brightness temperature and on the IC emission of the source 
were tested simultaneously against each other. 

The first results of the ground-based observing campaign 
are presented in Sect.|2lof this paper; the INTEGRAL observa- 
tions are described in Sect.|3j the simultaneous spectral energy 
distribution of the source is presented in Sect.|^ in Sect. [S] we 
discuss our results and draw conclusions. 



2. Ground-based multifrequency observations 

The multiwavelength observations, organized under the frame- 
work of the ENIGMA' collaboration together with a WEBT^ 
campaign, were scheduled for the period November 06-20, 
2003 (hereafter referred to as the core campaign) to provide 
the low-energy counterpart of the INTEGRAL observation of 
S5 0716+71. Since the source underwent an unprecedented 
outburst phase at radio and millimeter frequencies during 
September-October 2003, the campaign was started earlier 
(October 2003). When the core campaign began on November 
06 (JD = 2452949.5), S5 0716+71 was still brighter in the ra- 
dio and millimeter domain than in all past studies. A rather 
faint optical state (R ~ 14.17 - 13.64) characterized the source 
in this band for the whole duration of the INTEGRAL point- 
ing. During the last two days of the core campaign the source 
entered a brightening phase, and the ground-based monitoring 
was continued till May 2004. 

We here present the results of the core campaign in terms 
of radio light curves at 32 and 37 GHz, and an optical 7? -band 
light curve. Flux variation ranges at other radio, millimeter and 
sub-millimeter frequencies are included in the simultaneous 
SED presented in Sect. 0] the corresponding light curves will 
be discussed in Agudo et al. (2006 l and Fuhrmann et al. (in 
prep.). The ground-based observing facilities whose contribu- 
tion to the core campaign is presented in this paper are listed in 
TablelU 

2.1. Radio data 

2.1 .1 . Observations, data reduction and calibration 

The radio observations of S5 0716+71 at 32 GHz were per- 
formed during JD = 2452954.889-2452961.251 with the 100- 
m radio telescope in Effelsberg. The measurements were car- 
ried out using repeated scans, each of them consisting of 4-8 
cross-scans in azimuth and elevation, resulting in typical single 
flux-density exposure times of 120^80 s. The antenna tem- 
perature on the source was measured through the averaged 
pointing-corrected amplitude of the Gaussian-shaped cross- 
scans and with the standard data reduction software of the 100- 
m telescope. Regular measurements of the system temperature 
were used to determine the time-dependent fluctuation of the 
atmospheric opacity, which was used to correct the measured 
antenna temperature for each scan. Several non-variable se- 
condary calibrators (within 10°-20° of S5 0716+71) were ob- 
served with the same sampling of the target source, in order 
to ensure an accurate calibration of the residual gain fluctua- 
tions, mainly due to atmospheric effects. The total flux-density 
scale was fixed using the standard primary calibrators (e.g. 3C 
286, 3C 295, NGC 7027), whose fluxes are given in Baars et 
al. J 19771 1 and Ott et al. ( 1994 1. Further details on the observing 
method and data-reduction scheme can be found in Kraus et al. 
( 2003 ) and references therein. 

' http : //www . Isw . uni -heidelberg . de/pro jects/enigma/| 

- http : //www . to . astro . it/blazars/webt/] 
(see e.g. Mattox et al. il998. . Villata et al. .2004. Raiteri et al. l2005l and 
references therein) 
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Table 1. List of the ground-based observing facilities whose contribution to the core campaign is presented in this paper in the 
form of light curves and/or flux variation ranges in the SED. a: Radio, millimeter and submillimeter antennas, listed in ascending 
order of the lowest observing frequency, b: Optical telescopes, listed in order of longitude, and their contribution to the core 
campaign in terms of /?-band useful data (and observing nights). The offsets indicate the corrections applied to the data sets in 
order to eliminate instrumental inconsistencies (when a range is given, different offsets were used on different nights). 



a. Radio, millimeter and sub-millimeter observatories 


Observatory 


Location 


Telescope diameter 


Observing frequencies 


(telescope) 




(m) 




(GHz) 


Westerbork 


The Netherlands 


14x25 


1.392, 1.67 


UMRAO 


Michigan, USA 


26 


4.S 


i, 8.0, 14.5 


EfFelsberg 


Germany 


100 


4.85, 10.45, 32.0 


Metsahovi 


Finland 


13.7 


22.2, 36.8 


IRAM 


Spain 


30 


86.0, 230.0 


HHT-SMTO 


Arizona, USA 


10 




345.0 


JAC (JCMT) 


Hawaii 


15 


350.0, 664.0 


b. Optical observatories 


Observatory 


Location 


Telescope diameter 


i?-band data Offsets 






(cm) 


(nights) 




Lulin 


Taiwan 


100 


153 (7) 


0. 


Mt. Maidanak 


Uzbekistan 


150 


2(1) 


0. 


Abastumani 


Georgia 


70 


655 (6) 


[-0.03;0.] 


Crimean 


Ukraine 


70 


3(2) 


+0.03 


Tuorla 


Finland 


103 


102 (2) 


[+0.01;+0.02] 


MonteBoo 


Czech Republic 


62 


59 (2) 


[+0.05;+0.07] 


Perugia 


Italy 


40 


7(4) 


0. 


Heidelberg 


Germany 


70 


333 (8) 


[+0.03;+0.06] 


Michael Adrian 


Germany 


120 


90(4) 


+0.03 


Torino 


Italy 


105 


19 (3) 


-0.02 


Hoher List 


Germany 


106 


41 (2) 


+0.14 


Calar Alto 


Spain 


220 


374 (3) 


[+0.01;+0.02] 


Roque de los Muchachos (KVA) 


Spain 


35 


58(5) 


-0.01 


Roque de los Muchachos (WHT) 


Spain 


420 


751 (1) 


0. 


Bell 


Kentucky, USA 


60 


3(1) 


0. 


St. Louis 


Missouri, USA 


36 


3(1) 


0. 


Kitt Peak (WlYN) 


Arizona, USA 


90 


21(6) 


0. 


Coyote Hill 


California, USA 


28 


73(1) 


0. 


University of Victoria 


Canada 


50 


102(1) 


0. 



The measurements at 37 GHz were obtained in the period 
JD = 2452951.878-2452962.838 with the 13.7-m diameter 
radome-enclosed antenna of the Metsahovi Radio Observatory. 
The 37-GHz receiver is a dual horn, Dicke-switched receiver 
with a HEMT preamplifier, and it is operated at room tem- 
perature. On-on observations were performed, alternating the 
source and the sky in each feed horn, and adopting typical in- 
tegration times of 1200-1400 s. The source DR 21 was used as 
the primary flux calibrator, whereas 3C 84 was used as a secon- 
dary calibrator. Errors in the calibrated fluxes were computed 
by taking both the contribution of the measurement rms and the 
uncertainty in the absolute calibration into account. More de- 
tails about the Metsahovi observing system and data reduction 
can be found in Terasranta et al. ( I1998> and references therein. 

2.1 .2. The light curves 

FigH shows a superposition of the 32 GHz (black symbols) 
and 37 GHz (grey symbols) radio light curves. In spite of the 
different measurement accuracy and scatter which characterize 



the two data sets, a rise of the source flux over the observing 
period, with some shorter-term fluctuations superimposed, is 
clearly visible at both wavelengths. 

A variability test was performed on the two data sets by 
checking their consistency with a constant flux level at the 0.1% 
significance level: the test yielded a negative result, meaning 
that the light curves clearly exhibit variability. 

The overall increase of the brightness was approximated by 
a linear slope: the 32-GHz data can be represented by a linear 
increase of the flux of 41% in 6.4 days, whereas the 37-GHz 
data can be described by a flux rise of 39% in 10.7 days (see 
Table|2li). 

A consistency check of the two light curves was per- 
formed for the period in which they overlap, namely JD - 
2452954.878-2452958.962 . In this period, the data at 32 and 
37 GHz were characterized by a sampling of ~ 12 and ~ 15 
data-points per day, and average flux uncertainties of 4.6% and 
5.0%, respectively. The difference in scatter of the two data 
trains is reflected in the fractional root-mean-square variabil- 
ity amplitude Fvar (see Vaughan et al. 120031 and references 
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Fig. 1. Radio-optical light curve of S5 0716+71 during the core campaign. Panel a: 32 GHz radio light curve, and 37 GHz radio 
light curve scaled to the 32 GHz one (scaling factor: < /^32GHz/^37GHz >=0.89). Panel b: /?-band optical light curve. The shaded 
strip indicates the period of the INTEGRAL pointing. The two radio light curves, consistent with each other, exhibit significant 
inter-day variability but low-amplitude IDV. The optical light curve displays a weekly modulation of the emission with typical 
IDV features superimposed. 



therein), which is 9.1% at 32 GHz and 13.9% at 37 GHz. 
However, a Kolmogorov-Smirnov test showed that the diffe- 
rences between the 32-GHz and the 37-GHz data (rescaled to 
the 32-GHz frequency) are normally distributed at 95% signi- 
ficance level, confirming that the measurement errors can well 
account for the discrepancies between the two curves. The light 
curves are hence consistent with each other. 

The period of overlap of the two light curves covers the 
inter-day flux variation recorded in the source during the core 
campaign: the rise in flux is consistent at both frequencies with 
the linear increase of 35% in 4.1 days which best fits the 32- 
GHz data (see TablelS?). 

A search for faster variations was then performed in both 
light curves for the whole core-campaign period, during which 
the sampling of each light curve is such that it is easy to identify 
subsets of data separated in time by gaps of ~ 8 hours. 

At 32 GHz, four subsets of light curve can be fitted by a 
linear trend, whose slope is significantly variable. The result of 
the variability test (made by assuming the slope equal to zero) 



and the best-fit value of the slope are reported, for each sub- 
set, in Table Seven subsets of data at 37 GHz can also be 
approximated by linear slopes (see Table |2j;), with the excep- 
tion of seven measurements during JD = 2452952, 2452954, 
and 2452956 (see Table IJJl). The most extreme deviation was 
the one recorded during JD - 2452956 (which yielded a flux 
difference of ~ 4cr in the distribution discussed above), cor- 
responding to a linear increase of the flux of 42% in 0.12 days. 
The other two strong deviations can be described by linear flux 
variations of 33% in 0.41 days and 20% in 0.12 days, respec- 
tively. 

The reliability of the observations can be best estimated by 
comparing the measurements taken simultaneously with diffe- 
rent telescopes (see e.g. Witzel et al., 1992 for a comparison 
of the 5-GHz observations of S5 0716+71 performed at the 
same time with the VLA and the Effelsberg telescope in May 
1989). Therefore, although a careful analysis of the 37-GHz 
data showed that the measurements characterized by extreme 
deviations were not affected by any obvious systematic effect. 
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Table 2. Results of the linear fits to the radio Ught curves at 32 
and 37 GHz over different periods during the core-campaign 
(see text for details). Col 1: time interval relevant to the fit; 
Col. 2: observing frequency [and size of the time bin applied to 
the light curve in order to make the linear fit acceptable at 5% 
significance level]; Col. 3: linear slope; Col. 4: reduced of 
the fit; Col. 5: Fit significance. 



Time interval v (GHz) Slope x,-ed Sign. 



(JD-2452000.00) [bin (h)] (Jy/d) (%) 



Q^A QSQ Qfi^ 0^ 1 


10 

dL 


a 


_L.n 1 Q -1- n ni 


1 nn 


Al '^.l 
4 / .J J 


Q-^ 1 fi7Q Qfil Q^S 






+U. iz ± yj.Kjj 


1 ISl 

1 .JO 


1 Q C\A 


CiZA QHQ CiZQ O^^O 


JZ 


D 


+U.ZJ ± U.UZ 




JJ.JU 




3/ L3.J 




+U.OJ ± U.l J 


1 Tl 

l.zz 


n 1 
Zz. / I 


c\ZA oon nc^ Qm 


3z 


c 


n 
U. 


l.zi 


1 n m 




JZ 




1 n no 1 n n/i 


1 IT 
1.1/ 


T /I no 


nc/i riA/i nc7 n/zo 

950.91)4-957.968 


32 




0. 


0.86 


61.50 




JZ 




-1-0 -1- n 1 ^ 

-rU. J J It U. 1 J 


U.JO 


70. JU 


958.372-959.711 


32 




0. 


0.36 


99.67 




j2 




—U.UI ± \j.\jy 


U.Jo 


QQ A ^ 


you.oi^H — Vol .zj 1 






n 
U. 


U.oU 


/ 1 .JO 




32 




+0 48 + 28 


0.28 
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-0.10 + 0.37 
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1.16 




37 




0.05+0.17 


2.24 
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37 
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2.23 
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37 
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2.23 


6.30 
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37 
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+2.61+0.33 
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37 




-2.14 + 0.30 


2.32 


2.32 


956.225-956.428 


37 




+8.86+1.21 


2.88 


3.47 



the lack of any counterpart at other radio frequencies did not 
allow us to consider the lower limits to the brightness tempera- 
tures inferred from those variations (up to Tb ~ 1.5 x 10'^ K) 
as a reliable constraint on the variability brightness temperature 
of the source. 

2.2. Optical data 
2.2.1 . Observations 

Table^) (cols. 1^) shows names and locations of the nineteen 
observatories which contributed to the core campaign with R- 
band optical data, together with the sizes of the telescopes used 
and the amount of useful observations provided. 



It was recommended that the optical observers performed 
intranight observations of S5 0716+71 and of the field refe- 
rence stars nos. 1 to 8 whose magnitudes are given in Villata 
et al. ( 1998). The choice of the exposure times was left to the 
observers, who found the optimal compromise between high 
accuracy and good temporal sampling. 

The core-campaign period was characterized by un- 
favourable weather conditions in most of the observing sites. 
However, the high declination of the source, the long northern 
nights, and the large number of telescopes involved allowed a 
considerable overlap among observations carried out with tele- 
scopes located at very different longitudes. This led to unprece- 
dented dense monitoring of the source over such a long period: 
2849 useful observations were performed in ~ 15 days, with an 
average rate of ~ 8 frames per hour. 

2.2.2. Data reduction, calibration and assembly 

Ten out of nineteen data sets (44% of the collected data) 
were homogeneously analysed with aperture photometry on 
de-biased and flat-fielded frames. The remaining data sets were 
analysed with different procedures by the observers them- 
selves. 

The instrumental magnitudes were then processed in order 
to obtain the standard magnitudes of S5 0716+71 and the re- 
levant errors. Due to the variety of characteristics of the tele- 
scopes and detectors involved, as well as of the collected pho- 
tometric results, not all the reference stars suggested to the ob- 
servers were available for the calibration. Moreover, because 
the source was faint during most of the core campaign, the 
long exposure times required often caused the saturation of re- 
ference stars 1 and 2. In order to achieve a homogeneous data 
calibration and minimize the instrumental offsets, we defined 
our calibration sequence of reference stars common to most 
of the data sets and free of saturation, i.e. stars 3, 5 and 6. 
The uncertainty in the calibrated magnitudes of S5 0716+71 
was estimated, for each subset of data taken during the same 
night and under comparable observing conditions, as follows: 
the two stars of the calibration sequence closest in magnitude to 
the source were selected, and their magnitude difference calcu- 
lated. For each frame, the deviation of this quantity with respect 
to its mean was found. The error was taken as the larger of this 
deviation and the standard deviation of the difference over the 
night. 

Binning (over time intervals from 5 to 20 minutes) was ap- 
plied to intranight data trains affected by noise and/or low ac- 
curacy due to either non-optimal sky conditions or too short 
integration times. 

The calibrated light curves relevant to the various tele- 
scopes were then assembled according to the procedure de- 
scribed in Villata et al. (^002). Instrumental offsets were com- 
puted, night by night, whenever this was allowed by the tem- 
poral overlap between two data sets, i. e. when one of the two 
datasets had at least two data points belonging to the time inter- 
val covered by the second data set. Each offset was defined as 
the mean magnitude difference of the overlapping parts of the 
two light curves, the difference being computed between pairs 
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of data points separated in time by no more than 5 minutes 
(see Table[n5, col. 5). The data sets were then corrected by the 
above offsets. Data subsets characterized by both lower accu- 
racy and worse temporal sampling were finally discarded when 
higher-accuracy and better-sampled data trains were available 
in the same time interval. Data points affected by errors greater 
than 0.05 mag were not included in the final light curve, unless 
no other data were available within 10 minutes. The average 
accuracy of the final light curve is of order ~ 1%. 

2.2.3. The light curve 

The /?-band core-campaign light curve is displayed in Fig.^. 
Flux densities were derived from magnitudes using the absolute 
calibrations of Bessel (T979V and dereddened with the extin- 
ction laws of Cardelli et al. ( 1989 1, under the assumption of a 
Galactic extinction Ab - 0.132 mag (provided by NED; from 
Schlegel et al. [T998l . 

This light curve is characterized by a mean sampling of 20 
minutes through the whole core-campaign period, and of about 
14 minutes during the INTEGRAL observation. 

The source displayed remarkable variability during the 
core campaign: the overall peak-to-peak /?-band variation was 
AFr = 12.1 mJy (AR = 1.07 mag), and the fractional root 
mean square variability amplitude Fyar = 23.3%. 

A weekly modulation of the emission with intra-day fea- 
tures superimposed characterized the relatively faint state of the 
source before JD = 2452963.0: during this period the source 
brightness varied between 7.2 and 12.4 mJy (R - 14.17 - 
13.58), the mean flux was < Fr >= 9.5 mJy (< R >= 13.87), 
and Fvar = 10.6%. After JD = 2452963.0, a brightness rise 
of 7.8 mJy (AR - 0.57 mag) occuiTed in 25.7 hours, and the 
source reached 19.2 mJy [R - 13.10 mag), the highest level 
recorded during the core-campaign period. 

The IDV features do not generally show well-defined 
shapes, suggesting that there might be blending between suc- 
cessive flares. Their linear parts (in magnitude scale), lasting 
up to ~ 3.5 hours, display rising and declining rates not faster 
than 13% per hour These variation rates are comparable with 
the steepest slopes reported by Wagner et al. ( 1996 1 (~ 10% 
h~') and also with the slopes found by Villata et al. (I2000> 
(12% h"'). Lower-amplitude intra-day variations, on compara- 
ble time scales, were reported by Ghisellini et al. (I1997> . Sagar 
et al. {19991 and Nesci et al. <Q002\. Detailed temporal and 
spectral analysis of the optical variability will be given in forth- 
coming papers. 

3. INTEGRAL observations and data analysis 

INTEGRAL (Winkler et al. 2003 1 observed S5 0716-1-71 from 
2003 November 10th, 11:20:04 UT to 2003 November 17th, 
09:09:31 UT, for a total amount of ~ 539 ks. 

The observations were strongly affected by the biggest ever 
recorded solar flare (classified as X28), which occurred on 
November 04, 2003. Because of this event, SPI (Vedrenne et al. 
[2003 1 underwent annealing treatment during revolutions 132- 
136 and was not operational for most of the granted obser- 
ving time. Of the two detectors of JEM-X (Lund et al. .2003j . 



only JEM-X 2 was used, while JEM-X 1 was switched off. 
IBIS/ISGRI (Ubertini et al . (20041 Lebrun et al. l2003t and 
IBIS/PICsIT (Ubertini et al. l2004l Di Cocco et al. l2003l were 
both operational. 

The INTEGRAL data were analysed by means of the 
INTEGRAL Offline Scientific Analysis (OSA) software, 
whose algorithms are described in Westergaard et al. (I2003> 
for JEM-X and Goldwurm et al. (2003 1 for IBIS. Version 4.2 
of the OSA software was used for the analysis of the IBIS data, 
whereas the improved OSA 5.0 was necessary for the JEM-X 
data analysis. The PICsIT data could not be analysed due to the 
failure of the OSA software pipeline. 

3.1. JEM-X 

The JEM-X observation is split into 145 different science win- 
dows with an average exposure of ~ 3 ks each. The OSA soft- 
ware was used to look for sources within an offset angle of 5° 
to avoid spurious detections. S5 0716+71 was not detected in 
any of the single science windows, neither in the 3-35 keV band 
nor in any sub-band. There was also no other detected source 
in the field of view. 

We combined the images from all the individual science 
windows into a mosaic (coiTesponding to a total effective expo- 
sure of ~ 381 ksec) with the varmosaic routine of the FTOOL 
package; however, the source remained undetected. 

We estimated an upper limit to the flux of S5 0716+71 by 
means of the statistical variance of the mosaic intensity at the 
position of the source. The 3cr upper-limit intensity (cts cm 
s"'), with cr being the square root of the variance, was multi- 
plied by the ratio of the flux of the Crab Nebula in the JEM- 
X band (Toor & Seward 1974, N0rgaard et al. 1994) and the 
JEM-X intensity obtained from a mosaic image of the Crab. 
Assuming the spectrum of S5 0716+71 to be identical to that 
of the Crab (F = 2.1 + 0.03; see Toor & Seward 1974), we ob- 
tained a 3cr flux upper limit of ^3-35 keV = 6.12 x 10"'^ erg 
cm"^ s"' (see Table|3t), which is lower than the upper limit de- 
rived by Plan et al. (20Q5j from a shorter (189 ks) observation 
performed in April 2004. 

3.2. IBIS/ISGRI 

The ISGRI-instrument data set consists of 150 science win- 
dows, coiTesponding to a total effective exposure of 362 ksec. 
The standard energy binning was used for the image deconvo- 
lution, and the individual frames were accumulated into a co- 
added mosaic image. Several sources, including Mrk 3 and Mrk 
6, were reported as significant detections by the OSA software. 

S5 0716+71 was not detected by the source locator at 3o" 
significance in any energy bin of the mosaic image. A signal 
was detected at a significance greater than 1 cr in at least one e- 
nergy bin between 15 and 200 keV in 54% of the individual sci- 
ence windows. The reliability of these detections was checked 
by extracting, from each science window, the count rate at the 
sky position of S5 0716+71 and in 12 background regions 
located within 1° of the source. A two-sample Kolmogorov- 
Smirnov test showed that the resulting distributions of source 
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Table 3. Results of the INTEGRAL analysis: JEM-X and 
IBIS/ISGRI 3cr upper limits to the flux of S5 0716+71 during 
November 10-17, 2003. 



a. JEM-X 



Energy 
range (keV) 


Intensity 

(10"^ cts/cm-/s) 


Photon 
index T * 


Flux 
(10"'- erg/cmVs) 


3-35 


< 5.89 


2.1 


< 6.12 


b. IBIS/ISGRI 


Energy 
range (keV) 


Count rate** 
(cts/s) 


Photon 
index T * 


Flux 
(10"" erg/cm^s) 


15-40 
40 - 100 
100 - 200 


< 0.176 

< 0.160 

< 0.140 


1.6-2.1 
1.6-2.1 
1.6-2.1 


< 1.34- < 1.41 

< 1.70- < 1.72 

< 5.58- < 5.73 



* Photon index adopted to derive the flux from the count rate, under 

the assumption that F{E) oc E'^'^ 
** Count rate errors are Icr uncertainties 



and background count rates in each sub-band were consistent 
with them being drawn from the same parent distribution at 
the 5% significance level, confirming the non-detection of S5 
0716+71 at the science-window level also. 

In order to estimate an upper limit to the source flux in the 
ISGRI energy band, we extracted, for each sub-band, the value 
of the statistical variance of the count rate at the position of 
the source from the mosaic^ of all the science windows. The 
coiTesponding standard deviation (cTstat) was multiplied by the 
HWHM of the distribution of the significances of the mosaic 
(cTsyst) to take systematic errors into account (A. Goldwurm, 
priv. comm.). The 3cr upper limits (cr = CTstat ■ o"syst) to the 
count rate derived with this procedure are shown in Table (Sj). 

The above-mentioned result was crosschecked by estima- 
ting, for each energy sub-band, the standard deviation iXstat 
from the statistical variances of all the science windows (A. 
Goldwurm, priv. comm.). The flux and variance at the nominal 
position of the source were extracted from each science win- 
dow. The mean of these fluxes was then computed by weigh- 
ting it by the inverse of the corresponding variances. The error 
in the mean flux, representing iXstat, was found to have the same 
value as in the case of the mosaic, thus yielding equivalent up- 
per limits. 

The count-rate upper limits derived with these methods 
were converted into fluxes through XSPEC by using the re- 
sults of the spectral analysis of three BeppoSAX observations 
(Giommi et al. [T9991 Taghaferri et al. i2003) above ~ 2 keV, 
where the source spectrum can be well represented by a power 
law (T - 1.6 - 1.96). In addition, the Crab Nebula spectrum 
(r - 2.1) was also considered, for consistency with the as- 
sumption made for the computation of the JEM-X flux. By 

^ Note that the parameter responsible for pixel spread in con- 
structing the mosaic was not activated in our ISGRI analysis 
(OBSl J'ixSpread=0), in order to obtain a better estimate of variance 
in the mosaic (A. Goldwurm, priv. comm.) 



varying the assumed photon index in the range Y - 1.6-2.1, 
we obtained the 3cr flux upper limits given in Tabled?. 

A comparison of this result with the ISGRI flux of 
S5 0716+71 derived by Pian et al. (2005 1 in the 30-60 keV 
range from their ~ 30% shorter observation (256 ksec) of April 
2004 was performed by computing the flux in the same ener- 
gy bin, under the assumption of a spectral index equal to that 
of the Crab Nebula. The upper limit obtained in this way was 
^^30-60 keV = 1.06 X 10"" erg cm"^ s"', a factor of ~ 3 below 
their detection. 

4. The spectral energy distribution (SED) 

The broadband monitoring carried out during the INTEGRAL 
pointing (see Table [l] for a list of the ground-based observing 
facilities) allowed us to assemble, for the first time, the SED of 
S5 0716+71 with truly simultaneous data spanning more than 
10 decades of frequency and characterized by an exceptional 
energy sampling. 

The SED of S5 0716+71 is displayed in Fig. |2 whereas 
the Fy vs V diagram is shown in Fig. |3] (Online Material). The 
flux variation ranges at radio-to-optical frequencies recorded 
during the INTEGRAL pointing period (see Sect.|2}, as well as 
the upper limits estimated from the non-detection of the source 
by INTEGRAL (see Sect. |3} are plotted with dark diamonds, 
and are superposed on the historical data or variation ranges, 
drawn with grey symbols^. 

The low-energy part of the SED, which is commonly in- 
terpreted as due to synchrotron emission, shows the unprece- 
dented brightness of the source in the radio-millimeter ener- 
gy range and the moderate emission recorded at optical fre- 
quencies. The synchrotron peak (representing the most relevant 
output of the synchrotron component) is likely to be located 
around 10'^ - 10"* Hz. 

During the core campaign, the synchrotron spectrum of S5 
0716+71 was inverted {a > 0, with Fy oc V) in the radio- 
millimeter wavelength domain, with a turnover frequency close 
to 90 GHz (see Agudo et al. I2006t . The spectral index of 
the partially opaque part of the spectrum was derived from 
the simultaneous radio measurements at 5 and 32 GHz (con- 
siderably below the turnover) performed on JD=2452955 and 
JD=2452961 with the Effelsberg 100-m radio telescope (see 
Fuhrmann et al., in prep.) Its value is 05-^2 - +0.3 on 
JD=2452955, and 0-5-32 = +0.5 on JD=2452961: as flie 
millimeter-flux rose towards the end of the observing period, 
the radio spectrum became more inverted. 

As usually found in compact radio sources, the spectrum 
at frequencies lower than the turnover was not rising as fast as 
the optically thick spectrum of a homogeneous spherical syn- 
chrotron source (Fy oc v^'^). Such behaviour is traditionally in- 
terpreted as the result of a superposition of spectra characte- 
rized by diff'erent synchrotron self-absorption (SSA) frequen- 
cies, which might be produced either by a finite number of ho- 
mogeneous components (Kellermann & Paulinv-Toth I 1 9691 or 

* Previous studies of the SED of S5 0716+71 can be found in 
Wagner et al. jT9951 , Ghisellini et al. HWTli, Giommi et al. kl999t . 
Ostorero et al. 1200 1> . and Tagliaferri et al. 12003 J . 



L. Ostorero et al.: Testing the Inverse-Compton catastrophe scenario in S5 0716+71 



9 




10 15 20 25 

log V [Hz] 



Fig. 2. Spectral energy distribution (SED) of S5 0716+71: dark diamonds represent data simultaneous with the INTEGRAL 
pointing, and grey symbols represent historical data; variation ranges are indicated by vertical bars. Simultaneous data are from 
this work: they show the exceptionally bright state of the source in the radio-to-submillimeter domain, the moderate level of the 
optical emission and the upper limits to the hard-X-ray brightness. Historical data are from: RATAN-600, at the wavelengths of 
1.38,2.7, 3.9,7.7, 13 and 31 cm (Ostorero et al., in prep.); Kuhr et al. ( 1981 1, Waltman et al. (1981 1, Eckartet al. (1982i, Perley 
l lT982t . Perley et al. (fT982t . Lawrence et al. ( fT985t . Saikia et al. ( 1987|, Kiihr & Schmidt ( 1990 1, Moshir et al. ( 1990|, Hales 
et al. i tTWTt . Krichbaum et al. (H^^, Gear et al. (TW4i, Hales et al. fW^, Douglas et al. ( 1996 1, RengeHnk et al. ( ITWTt . 
Zhang et al. (jl997 1, Riley et al. ( 1999l, Cohen et al. (2002 1, Raiteri et al. (2003 ) and references therein, at other radio-to-optical 
frequencies; Pian & Treves ( 1993 1, & Ghisellini et al. ( 1997 1 in the UV band; Biermann et al. | |1992> . Comastri et al. ( I1997> . 
Kubo et al. dT^^ . Giommi et al. ( 1999 1, Tagliaferri et al. (l2003t . & Pian et al. (I2005i) in the X-ray band; McNaron-Brown et al. 
| |1995> . Hartman et al. Cl999j . & Collmar (■2006j in the y-ray energy range. 



by an inhomogeneous source with gradients in the magnetic 
field and particle density (Condon & Dressel [19731 de Bruyn 
[T976l Marscher[T977i. 

The variability observed in the radio-millimeter energy 
range, of which we showed an example in the 32 and 37 GHz 
light curves presented in Sect. 12.11 hence occurred in a regime 
in which SSA processes might have played an important role. 

The IDV study carried out during the fainter radio-to- 
optical state of the source of February 1990 (Wagner et al. 
fl990 Quirrenbach et al. ll991l l. reported a variable radio spec- 
trum which was on average optically thin both in the 5.0-8.4 
GHz energy range (cs-s^ - -0.10) and in the 1.4-5.0 GHz 
band (ai 4-5 = -0.35) (Wagner et al. I1996t . The long-term 
behaviour of the radio spectral indices indicates spectra on 
average inverted in the 4.8-14.5 GHz band (Q'4.g-8 - 0.72; 
org-u.s = 0.16) during 1985-1992 (Wagner et al.|T5^, as well 
as in the 5.0-15.0 GHz band (ag-is = 0.19) during 1978-2002 
(Raiteri et al. 20031, where also a flatter-when-brighter trend 
and a weak correlation between spectral flattening and ejection 
of a new component were recognized by Raiteri et al. (2003 1 
and Bach et al. (■2005j . respectively. However, no other truly 



simultaneous spectral information on the source spectrum from 
radio up to millimeter frequencies is available to date. 



The moderate level of the optical emission recorded during 
the bright radio state is consistent with the absence of correla- 
tion, at zero time lag, between major optical and radio events in 
the source, as already noticed by Raiteri et al. (2003) from in- 
spection of the historical ( 1 978-2002) radio-optical light curve. 



As far as the high-energy part of the SED is concerned, 
the upper limits provided by JEM-X (see Sect. 13. 1> and by 
IBIS/ISGRI (see Sect. 13. 2> in the two lower-energy bins are 
consistent with the levels recorded by ASCA (Kubo et al. 
0998) and BeppoSAX (Giommi et al.n^51 Taghaferri et al. 
|2003 ) during past observations, and indicate that the source was 
fainter than during the INTEGRAL observation of April 2004 
(Pian et al. 2005 1, represented in Fig.|2lby the grey triangle at 
~ 41.4 keV (~ 10''' Hz). The two higher-energy IBIS/ISGRI 
upper limits are comparable with the limits derived from the 
OSSE observations (McNaron-Brown et al. ri995) . 
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5. Discussion and conclusions 

The known high duty cycle of S5 0716+71 was confirmed 
by the intensive simultaneous radio-to-optical monitoring per- 
formed during the core campaign. The source was variable in 
both the radio and optical bands, although the amplitude of IDV 
recorded in the radio regime was significantly lower than in 
previous studies. 

Let us assume long-term stationary radio emission, and 
consider a Gaussian brightness distribution for the variable 
source region. If the event which triggers the flux fluctuations 
propagates isotropically through the source, the flux variability 
observed in the radio band allows one to estimate the brightness 
temperature of the source through the relation: 



Tb = 8.47X lO'^Fv 



A dh 



(1) 



(3) 



where Fy is the flux density in Jy, A is the wavelength in cm, ty 
in years is a function of the mean flux and the flux variation in 
the time interval considered according to the relation 

<Fy> ^ 

AFy (1+z) 

(Wagner et al.'1996', see also Jones et al. ll974> . z is the source 
redshift, and di^ is the luminosity distance of the source in Mpc, 
given by 

no Jo 

(Carroll et al. \T99% in the case of a A-dominated universe 
(Qa = 0.7, Q.M = 0.3, Qk = 0; see Spergel et al. 2003 1. 
Assuming Ho=72 km sec ' Mpc ' (Freedman 1200 U and a 
redshift z > 0.3, the source luminosity distance would be 
greater than 1510 Mpc. The brightness temperature derived 
from the overall increase of the source flux recorded at 32 
and 37 GHz during the period of overlap of the two-frequency 
measurements (At - 4.1 days) is hence determined to be 
Tb > (2.1 + 0.1) X 10"* K, if die observed flux variations 
are not affected by any propagation effects. Higher lower lim- 
its to the brightness temperature can be derived under less- 
conservative assumptions: a uniform distribution of the source 
brightness, and/or an anisotropic propagation of the pertur- 
bation which causes the variability would yield values up to 
Tb > (1.2 ± 0.1) X 10'^ K. Moreover, a redshift limit z > 0.52 
would increase Tb by a factor of ~ 3. In any case, our bright- 
ness temperature exceeds the IC limit of ~ 10'^ K by at least 
two orders of magnitude. 

If relativistic boosting of the radiation is the explanation of 
this excessive value, Doppler factors 6 <: (Tb/lO'^ K)'^-'(l+z) <: 
8 would be required in order to lower the intrinsic brightness 
temperature of the source below the theoretical limit^ . 

^ A more accurate value of the IC limit, and consequently a better 
estimate of the Doppler factor derived from the brightness tempera- 
ture, would require detailed knowledge of the characteristics of the in- 
trinsic spectrum of the source, including the synchrotron upper cutoff 
frequency, the self-absorption frequency, the radio spectral index, etc. 
(see e.g. Kellermann & Paulinv-Toth i 1 969i Blandford il990t Readhead 
fT994l. 



A recent analysis of the historical data set of VLBI images 
of S5 0716+71 at several frequencies performed by Bach et al. 
('2005 ) showed that Doppler factors ~ 20-30 are consistent with 
the proper motion of the jet components. These values, which 
rule out slower kinematic scenarios of the source jet presented 
previously (Witzel et al. 1988, Gabuzda et al. 1998), could 
largely account for the above brightness temperature, provided 
that the source region responsible for the observed flux varia- 
bility had the same kinematic properties as the VLBI jet com- 
ponents. 

The source flux evolution recorded in the centimeter- 
millimeter domain was always dominated by that of a core 
of size ^0.1 mas. Confirmation of this evidence was recently 
provided by the 15-GHz VLBA measurements during August 
2003 (Kovalev et al. 12005 1 and by the analysis of space- VLBI 
(VSOP) observations performed in 2000 (Bach et al. l2006t . 
The observed variability hence originates from either the sub- 
parsec jet or regions at the base of the jet itself. 

If the variability originates from one or more jet compo- 
nents moving according to the kinematics described by Bach 
et al. f'2005 1, and characterized by inverted or flat radio spec- 
tra (possibly becoming steeper as the components move away 
from the core), the brightness temperature would only appa- 
rently exceed the IC limit. The non-detection of the source by 
INTEGRAL might hence easily be the consequence of the non- 
occurrence of any Compton catastrophe in the source. This sce- 
nario would also enable us to reconcile the excessive VLBA 
brightness temperature of 7'b.j=o > 1.85 x lO'-' K derived by 
Kovalev et al. ( 2005 1 with a value lower than the IC limit, as- 
suming a Doppler factor 6 > (Fb/lO'^ K) (1 +z) > 24 (note that 
the VLBA brightness temperature scales by 5/(1 + z)). 

On the other hand, if the variable core emission is affected 
by a Doppler enhancement different from that of the resolved 
VLBI jet components, the IC-limit violation would be real 
for Doppler factors 6 ^ 8. In this case, Compton catastro- 
phes would have occurred in the source: the high-energy non- 
detection would hence provide a constraint for any model of 
the source emission taking second-order Compton scattering 
(e.g. Bloom & Marscher 1996 1, and hence the possibility of 
Compton catastrophes, into account. 

Some alternative explanations for the excessive brightness 
temperature take propagation effects into account. The corre- 
lation between optical brightness and radio spectral index, as 
well as the simultaneous change of variability time-scale ob- 
served during the February 1990 campaign (Wagner et al. 119901 
Quirrenbach et al. ll991l Wagner et al. 1996 1, have been among 
the strongest arguments against the extrinsic origin of radio 
IDV in S5 0716+71 in the past decade. However, the radio 
and optical intra-day variability observed during our campaign, 
when the source was in a brighter radio-to-optical state and had 
an inverted spectrum up to millimeter frequencies, do not ap- 
pear to be obviously coiTelated. This might be the consequence 
of optical-depth effects (which can modify the shape and am- 
plitude of flares) on radio emission intrinsically correlated with 
that in the optical band; alternatively, the radio and optical ra- 
diation observed might have come from non-cospatial compo- 
nents characterized by different sizes. In the absence of corre- 
lated radio-optical variations, the possibility of a contribution 
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of ISS to the observed radio variability cannot be completely 
ruled out, implying that the brightness temperature inferred 
from variability is not representative of the photon density of 
the source; in particular, the violation of the IC limit would 
again be only apparent in an ISS-driven variability scenario. 

It is very unlikely that this variability can be explained by 
scintillation alone. In fact, ISS is not a very efficient mechanism 
in the 32-37 GHz (8-9 mm) regime, to which our brightness 
temperature refers. Moreover, the simultaneous observations at 
86 GHz (3 mm), where ISS is less efficient than at 32-37 GHz, 
showed that ISS can be definitely ruled out as an explanation 
for the observed (comparable) variability at these frequencies 
and for the corresponding excessive brightness temperatures 
(Agudo et al. 2006i). Therefore, if ISS contributed to the ob- 
served 32-37 GHz flux evolution, in all likelihood it was not 
the dominant variability mechanism. 

In conclusion, violations of the brightness-temperature IC 
limit were inferred from the radio variability observed at 32- 
37 GHz. If the recorded radio emission is either intrinsic and 
beamed with Doppler factors 5 ^ 8, or strongly affected by 
ISS, the violation of the IC limit would be apparent, and readily 
consistent with the non-detection of the source by INTEGRAL 
in the X-y-ray regime. Intrinsic flux variations in the pres- 
ence of lesser beaming effects would instead imply a real vi- 
olation of the theoretical limit; in this case, the non-detection 
of corresponding y-ray avalanches provides a strong constraint 
for modelling of the Compton catastrophe in this source. At 
any rate, our unprecedented simultaneous broad-band measure- 
ments will help to define the parameter space of SSC emission 
models in detailed studies of the multifrequency properties of 
S5 0716-H71. 
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Fig. 3. Broadband spectrum of S5 07 16+7 1 : dark diamonds represent data simultaneous with the INTEGRAL pointing, and grey 
symbols represent historical data; variation ranges are indicated by vertical bars. See caption of Fig. 2 and Sect. 0] for more 
details. 



